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The article contains new results for spin-3/2 and -1/2 resonances. It specializes to the second 
resonance region, which includes the Pn(1440), Z?i3(1520) and 511(1535) resonances. New data on 
electroproduction enable us to determine the vector form factors accurately. Estimates for the axial 
couplings are obtained from decay rates of the resonances with the help of the partially conserved 
axial current (PCAC) hypothesis. We present cross sections to be compared with the running and 
future experiments. The article is self-contained and allows the reader to write simple programs for 
reproducing the cross sections and for obtaining additional differential cross sections. 
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I. INTRODUCTION 

In previous articles Q, @] we described the formal- 
ism for the excitation of the P33(1232) resonance. In 
the meanwhile, we extended the analysis to the sec- 
ond resonance region, which includes three isospin 1/2 
states: Pu(1440), P>i 3 (1520), Sii(1535). The dominant 
P33(1232) has been observed in neutrino reactions and 
there are several theoretical articles which describe it 
with dynamical calculations based on unitarized ampli- 
tudes through dispersion relations yL phenomenological 
El IE an d quark models B, M, llOj , as well as mod- 
els incorporating mesonic states including a cloud of 
pions. Articles in the past five years have taken a closer 
look at the data by analysing the dependence on neu- 
trino energy and Q 2 Q, 0, HJ . So far a consistent picture 
emerged to be tested in the new accurate experiments. 

For the higher resonances there are several articles, 
that describe their excitation by electrons 0, El El El, 
and only one by neutrinos. Experimental data for neu- 
trino excitation of these resonances are very scarce and 
come from old bubble-chamber experiments 0, El ITsl 
ITU |20 | . In the new experiments, studying neutrino os- 
cillations, there is a strong interest to go beyond the QE 
scattering 0, ^ and understand the excitation of these 
resonances. One reason comes from the long-baseline ex- 
periments where the two detectors (nearby and faraway) 
observe different regions of neutrino fluxes and kinematic 
regions of the produced particles. 

A basic problem with resonances deals with the deter- 
mination of their form factors (coupling constants and 
Q 2 dependences). The problem was apparent in the A 
resonance where after many years several of the form fac- 
tors and their Q 2 dependence became accurately known 
and were found to deviate from the dipoles. The sit- 
uation is more serious for the higher resonances where 
the results of specific models are used. In this article we 
adopt the approach of determining the vector couplings 
from helicity amplitudes of electroproduction data, which 
became recently available from the Jefferson Laboratory 
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[M l2J,|2a| and Mainz accelerator [2jj . This requires that 
we write amplitudes for electroproduction in terms of the 
electromagnetic form factors and then relate them to the 
vector form factors that we use in neutrino reactions. 
The above approach together with CVC uniquely speci- 
fies the couplings and the Q 2 -dependences in the region 
of Q 2 where data is available, that is for Q 2 < 3.5 GeV 2 . 

The axial form factors are more difficult to determine. 
For the axial form factors we adopt an effective La- 
grangian for the R — > Ntt couplings and calculate the de- 
cay widths. For each resonance we assume PCAC which 
gives us one relation and a second coupling is determined 
using the decay width of each resonance. 

Having determined the couplings for the four reso- 
nances, we are able to calculate differential and inte- 
grated cross sections. This way we investigate several 
properties in the excitation of the resonances. We find 
that a second resonance peak with an invariant mass be- 
tween 1.4 GeV and 1.6 GeV should be observable pro- 
vided that neutrino energy is larger than 2 — 3 GeV. Cal- 
culating cross sections in terms of the resonances provides 
a benchmark for their contribution and allows investiga- 
tors to decide, when more precise data becomes available, 
whether a smooth background contribution is required. 
Integrated cross section already suggest the presence of 
a background. 

In sections ITT1 IIIII and II VI we present the formalism and 
give expressions for the vector form factors. Estimates for 
the axial form factors are presented in the Appendix lAl 
Section El points out that the structure functions W4 
and W5 are important for reactions with tau neutrinos. 
We analyze differential and integrated cross sections in 
section IVTI We discuss there the existing data and point 
out a discrepancy in the normalization to be resolved in 
the next generation of experiments. 



II. ELECTROPRODUCTION VIA HELICITY 
AMPLITUDES 

One of the main contributions of this article is the de- 
termination of the vector form factors for weak processes. 
Our work relates form factors to electromagnetic helicity 
amplitudes, whose numerical values are available from 
the Jefferson Laboratory and the University of Mainz. 
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Values for the form factors at Q 2 = are presented in 
TableQ] Later on we express the weak structure functions 
in terms of form factors. 

In applying this method we must still define the nor- 
malization of electromagnetic amplitudes, which is done 
in this section. Before we address this topic we discuss 
the kinematics, the polarizations and spinors entering the 
problem. 

We shall calculate Feynman amplitudes shown in FigQ 
We estimate the amplitudes in the laboratory frame with 
the initial nucleon at rest and with the intermediate pho- 
ton moving along the z— axis. 




R,p'' 1 



FIG. 1: Electroproduction of the resonance. 

We define the four-momenta 

p" = (m*; 0,0,0), f = (<z°; 0, 0, q z ) 

and p ni =p t " + q" = (q° + m N ;0,0,q z ). 

The intermediate photon or W— boson can have three 
polarizations ej^ defined as 



-!'- m = +-^(0;l,±i,0), 



^ = -^0,0,,°) 



(III) 



and Q 2 = —q 2 . The spinors for the target nucleon, nor- 
malized as u(p, s z )u(p, s z ) — 2mjv, are given by 



i(p, s z ) = V2m N ( "q 2 ^ 



(II.2) 



where u, can be either 



U +l/2 



'-1/2 



The final resonances will have spin 1/2 or 3/2. For spin- 
3/2 resonances we shall use Rarita-Schwinger spinors 
constructed as the product of a polarization vector 
with a spinor U. The states with various helicities are 



defined by 

4 3/a) = 4 fl) w+V2) 

^ 1/2) = \flei S) U(p', +1/2) + ^ei R) U(p', -1/2) 
4~ 1/2) = \flei S) U(p>,-l/2) + y/iePutf, +1/2) 
4" 3/2) =ePu(j/, -1/2) 
with the spinor given as 



(II.3) 



U(p', Sz ) = ^Jp'° + M R 



and the polarization vectors by 



p' -3 



e?R,L) = +^(0;l,±i,0), efo = J-(p";0,0 )P '°). 

We emphasize that refer to the intermediate photon 
and belong to the J — 3/2 spin state. For spin- 1/2 
resonances the spinors are 



u{p',s z ) = ^p'°+M h 



U s 
p' -a 



p>°+M R Us ' 



(II.4) 



With this notation we can calculate three helicity am- 
plitudes for the electromagnetic process. For instance, 
for the D13 resonance the amplitude (R,+l/2\J em ■ 
£(r-)\N, — 1/2) in terms of form factors is presented in 
eqs. (jTOl and l|TVT6|l . 

We now define the overall normalization. Analyses 
of electroproduction data give numerical values for cross 
sections at the peak of each resonance 

sunups 



cjt(W = M r ) 



2lJlN ( A 2 + A2 \ 

M R Y R {A ^ 2 + A ^ 



a L (W = M R ) 



2m N Q 2 g2 

m r t r qi 



(II.5) 



(II.6) 



These are helicity cross sections for the absorption of 
the "virtual" photon by the nucleon to produce the final 
resonance mil. They are defined as 

a (l) {W) = lKj2\(RiH£(i)J?\N,s) 2 R(W,M R ) 



with K = 



4n 2 a 



W 2 -m 2 N 



The last factor in the cross section is the Breit-Wigner 
term of a resonance: 



R(W,M R ) 



MrTr 



1 



tt (W 2 -M 2 R ) 2 + M 2 T 2 R 



(11.7) 
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For a very narrow resonance or a stable particle it reduces 
to a 5— function. 

Numerical values have been reported for the ampli- 
tudes in (|II.5|1 . (|II.6(I which arc related to the following 
matrix elements 

(II.8) 

which we calculate in this article. 

Following standard rules for the calculation of the ex- 
pectation values the signs in these equations are deter- 
mined. There is an ambiguity for the sign of the square 
root, which we select for all resonances to be positive. 
Later on we must also select the sign for axial form fac- 
tors. We shall choose them in such a way that the struc- 
ture functions W3 for all resonances are positive, as in- 
dicated or suggested by the data. As a consequence the 
neutrino induced cross sections are larger than the cor- 
responding antineutrino cross sections. 



Comparing with l|IH.2|) . one easily sees, that, for 

the isospin-1/2 resonances, the weak amplitude satisfies 
the equality A(W+n -> = A\- A\. Since the 

amplitudes are linear functions of the form factors, the 
weak vector form factors are related in the same way to 
electromagnetic form factors for neutrons C™ and protons 
Cf: 

1=1/2: CY = C?-Cf, (III A) 

with index i distinguishing the Lorenz structure of the 
form factors. 

For the isospin-3/2 resonances one gets A 3 % (W + n — > 
i?(3)+) = A 3 p (W~p R^°) = y^a 3 . The weak form 
factors, which are conventionally specified for these two 
processes, are 

7 = 3/2: Cf = Cf = Cf. (111.5) 

For the process W + p — > R^ ++ the amplitude is \/3 
times bigger: A(W+p -> R ( ^++) = V3A(W+n 
i?( 3 )+). Some of the above relations were explicitly used 
in earlier articles 0,0,0. 



III. ISOSPIN RELATIONS BETWEEN 
ELECTROMAGNETIC AND WEAK VERTICES 

Our aim is to relate the electromagnetic to weak form 
factors using isotopic symmetry. The photon has two 
isospin components \I, I3) = |1,0) and |0,0). The isovec- 
tor component belongs to the same isomultiplet with the 
vector part of the weak current. Each of the amplitudes 
^3/2, ^1/27 "5*1/2 c & n be further decomposed into three 
isospin amplitudes. Let us use a general notation and 
denote by b the contribution from the isoscalar photon; 
similarly a 1 and a 3 denote contributions of isovector pho- 
ton to resonances with isospin 1/2 and 3/2, respectively. 
A general helicity amplitude on a proton (A p ) and neu- 
tron (A n ) target has the decomposition 



A p = A p (~f P -> 



A n 



A n ( in ^R°) = b + 



r 



|a 3 .(III.l) 



For the weak current we have only an isovector com- 
ponent of the vector current, therefore the b amplitude 
never occurs in weak interactions. A second peculiarity 
of the charged current is that V\ ± 1V2 does not have 
the normalization for the Clebsch-Gordon coefficients, it 
must be normalized as (Vi ± iV^/V^, which brings an 
additional factor of y/2 to each of the charged current in 
comparison with the Clebsch-Gordon coefficients: 

A(W + n^ R^+) = ^=a\ (III.2) 
v3 



IV. MATRIX ELEMENTS OF THE 
RESONANCE PRODUCTION, FORM FACTORS 

Following the notation of our earlier article 2], we 
write the cross section for resonance production in a no- 
tation close to that of DIS, that is we express it in terms 
of the hadronic structure functions. In the present nota- 
tion the sign in front of W3 in the cross section (lIV.lll is 
plus, 



da" 



N 



G 2 



W 



dQ 2 dW 4tt 



— COS t>c 



m 

m 2 N 



thnE 2 
1 



W^+m 2 ) 



2(k-p)(k'-p)-^m 2 N (Q 2 + ml) 



' N 



2 1 2 \ 



Wi , (Q 2 + ml) 



in 



2~ m l 

N 



-2^-m 2 Jk-p)\, (IV.l) 



■m 



N 



which implies that cross section for a reaction with neu- 
trino exceeds that with antineutrino if the structure func- 
tion W3 is positive. The corresponding formula for elec- 
troproduction is obtained by replacing the overall factor 
G 2 cos 2 Oc/^k by 2-na 2 /Q 4 and using the electromagnetic 
structure functions Wf 
In this case the contribution from Wf 
ligible and W| m = 



' m and yV| m instead of weak ones. 
em and Wl m 



is neg- 



A(W + p i? (3)++ ) = V2a 3 , 



A(W + n -> R 



(III.3) 
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TABLE I: Vector and axial couplings for the excitation of 
resonances (Q 2 = 0) 



R 


M R , GeV 


r tot , GeV 


elast 


[Jit NR. 


Cf 


cX 


CX 


cX 


P 33 (1232) 


1.232 


0.120 


0.995 


15.3 


1.2 


2.13 


-1.51 


0.48 


Ois(1520) 


1.520 


0.125 


0.5 


19.0 


-2.1 


-4.08 


1.51 


0.31 


R 


Mr, GeV 


Vtot, GeV 


elast 


g-rrNR 


9t 


aX 


gX 




Pn(1440) 


1.440 


0.350 


0.6 


10.9 


-0.51 


-4.6 


1.52 




Sii(1535) 


1.535 


0.150 


0.4 


1.12 


-0.21 


-4.0 


-1.68 





In the following subsections we specify the structure 
functions Wi for each resonance by relating them to form 
factors. To give the reader a quick overview of the results, 
we summarize the couplings (the values of the form factor 
at Q 2 = 0) in Table HJ One should keep in mind, how- 
ever, that all form factors have different Q 2 -dependences, 
which are given explicitly in the following subsections. 



current for the resonance production is expressed as 



(D 13 \r\N)=4 D \p>)d% i3 u(p) 



(IV.2) 



A. Resonance Z?i3(1520) 

We begin with a D13 resonance, which has spin-3/2 
and negative parity. The matrix element of the charged 



with u{p) the spinor of the target and the Rarita- 
Schwinger field for a -D13 resonance. The structure of dp 
is given in term of the weak form factors 



J 



cx_ 

m N 



in- 



cX 



cX 



—2-7 + — §-j> + — 5-p 

m N m z N m% 



m N 



III 



N 



75 - q 



C£_ 

m N 

r 



c£ 



'■N 



75 



q X q" 



C£ 



75- 



(IV.3) 



The general form of the current for D13 differs from 
that of P33 in the location of the 75 matrix, which is due 
to the parity of the resonance. The form factors Cf {Q 2 ) 
and Cf(Q 2 ) refer now to any D13 resonance. Later we'll 
specify them for the Di3(1520). The vector form factors 
are extracted from the clcctroproduction data, in partic- 
ular from the helicity amplitudes. We use recent data 
from l25l |30j, which were kindly provided to us by I. 
Aznauryan. 

Helicity amplitudes are expressed via the J em - e, which 
can be obtained from 



ljIVT2|) and (|IVT3|) by setting the 
axial couplings equal to zero and replacing the vector 
form factors by the electromagnetic form factors. This 
results in the following expression 



(D 13 \J em ■ e\N) = 4>^ m) F^u{N), 



with 



c 



(em) 



c 



(em) 



c, 



(em) 



m N 



'N 



N 



q"e h 



(IV.4) 



As it was discussed in the previous section, the elec- 
tromagnetic form factors of -D13 resonance are different 
for proton and neutron. Substituting the matrix element 
(|IV.4|I into Eqs. Ijll.8|l and carrying out the products 
with the spinors and Rarita-Schwinger field we obtain 
the following helicity amplitudes for electroproduction 



^3/2 



a 



(em) 



m N 



-(M R - mjv) 



a 



(em) 



III 



-q-p + 



c 



(em) 



-q-p 



N 



(IV.5) 
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.4 



D 13 
1/2 



C 



m 



3 (M H 



3 mpf 

(em) 

— q-p + 



2m 



a 



(ern) 



"""IS Hz \ 
m N 7T n ) 

M RP '° + M R ' 
-q-p (IV.6) 



A' 



'TV 



Dy = (1 + Q 2 /My) 2 denotes the dipole function with 
the vector mass parameter My — 0.84 GeV. To give an 
impression, how good this parametrisation is, we plot in 
Fig0thc hclicity amplitudes, obtained with these form 
factors. Vector form factors for the neutrino-nucleon in- 
teractions are calculated according to Eq. (|111.4fl 

For the axial form factors we derive in the Appendix 1X1 



qD 13 
J l/2 



2N q z 
3~Mfl" 



/i(em) n(em) 

63 (-M«) + 64 



m N 



in 



(Q 2 



2mjvt7° - m 2 N ) 



Ck 



where N 



200 
150 
100 
50 

-50 
-100 
-150 



m N 
JO 



-(q° 



N 



m N 



(IV.7) 



itin{W 2 — m N ) 



2m N (p n> + M R ). 



> 

o 



C\J 

<f 

C\J 
CO 

< 



A (P) 

A P 

A 1/2 

q (P) 
& 1/2 



* -i- 



0.5 1 



1.5 



a 



2.5 



3.5 



FIG. 2: Helicity amplitudes for the Di3(1520) resonance, cal- 
culated with the form factors from Eq. llV.8ll . 



Comparing Eqs. HIV. 51) . <|IV-6|) . I|IV.7I) for each value 
of Q 2 with the recent data on helicity amplitudes p3ll25l 
l3Cj . we extract the following form factors: 



Di 3 (1520) 



2.95/A 



V/8.9M V ' 



0) _ 



C 



-1.05/ D\ 



1 + Q 2 / 
-0.48 



.9M 2 ' 



4 n) = 

C< n) = 



-1.13/ Di 



(iv.; 



1 + Q 2 

0.46/Dy 



9M 2 '' 



l + Q 2 / 
-0.17 



.9M 2 



This is a simple algebraic solution with the numerical 
values for the form factors being unique. The function 



r A(D) 



m 2 N 



nl + Q 2 



C^ D \0) = -2.l (IV.9) 



Two other form factors and the Q 2 behaviour of the C§ 
can be determined either experimentally or in a specific 
theoretical model. To check how big the contribution of 
the C% and C± could be, we set them C% = C± = 1 
and computed in Fig(3] the various contributions to the 
differential cross section for E v = 2 GeV. Motivated by 
the results on P 33 (1232) resonance Q, the Q 2 dependence 
in our calculations is taken as 



C A(D) = 



C^ D \i))/D A 
1 + Q 2 /3M\ ' 



(IV.10) 



where Da = (1 + Q 2 /Ma) 2 denotes the dipole function 
with the axial mass parameter Ma = 1-05 GeV. 



2e-40 



> 
(3 



O 
-a 

~0 



1e-40 



5e-41 
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(c 4 A ) 2 
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:(c 5 A ) 2 \ 




\ c 3 A c 5 A 









0.5 1 

Q 2 , GeV 2 



1.5 



FIG. 3: Contribution of various form factors to the differential 
cross section. 

We conclude from Figure |31 that the contribution of 
(C4 1 ) 2 and C^C^ are small, but the other terms could 
be sizable. Their importance depends on the relative 
signs. It is possible, that C^Cf and C^C^ are positive 
and together with (C3 4 ) 2 give an additional 100% to the 
Z?i3 cross section. In case that C^Cf and/or C^C^ 
are negative, there are cancellations. In the following 
calculations of this article we set for simplicity C% = 

The hadronic structure functions for D13 resonance are 
similar to those for P33 , presented in our earlier paper Q 
and can be obtained from them by replacing m^M R by 
—rri]\[M R . We repeat the corresponding formulas here 
including the terms with and C3 4 , which could be 
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nonzero. The structure functions have a form where R(W, Mr) was defined in l|II.7|) and Vi are given 

2 below. In the following equations the upper sign corre- 

WdQ 2 ,v) = V l {Q 2 , v)R(W, Mr) (IV. 11) sponds to P 33 resonance and the lower sign to the L>i 3 . 



V l = ■ 1 1 ^[(Q 2 -q-p)Hq-P + rn%) + M R ttq-p) 2 + Q 2 m 2 N ±Q 2 rn N M R )] 

N R. 

[(Q 2 - q ■ P) 2 (q ■ P + m%) + M%((q ■ pf + Q 2 m 2 N T Q 2 m N M R )] 

jlVl D 

[(Q 2 -q-p)(q-p + m 2 N T 2m N M R ) - M R q ■ p] 



+ 



+ 



" l N li± R 

CYCY(Q 2 -q-p)- C\ C\ q • V 

m%M R 

C£Ct{Q 2 -q-p)- CiC^ [{q2 _ q . p){q . p + {Mr ± mjy)2) _ Af , q2] 



m 3 N M R 



[cX{Q 2 -q-p)-cU-p\ 



in 



\q-p + m 2 N =p m N M R ) + 



N 



a 



A Cf{Q 2 -q-p) 



(IV.12) 
[q- p + m 2 N ± m N M R ] 



V 2 



(Cr) 2 + (Cf) 2 
Ml 



cXc v 



Ar<A 



>2 \q-p + m 2 N + M 2 R ] + 3 J Q 2 [q.p+ {M R t m N ) 2 } + _ 3 J Q 2 [q ■ p + {Mr ± m N ) 21 



m N M R 



m N M R 



CXCX -Q 2 [q-p+{M RT m N ) 2 + Q 2 ] + + iC¥) y 2 J 2 M2R) + ^^SH ( - >2 \" ' " + '"a : "' ^ } 



m N M R 



+ ° 3 ° 5 m; q 



N 



1 %M R 



in 



N 



(r<A\2 TO jy , {Cj) 2 



Q 2 



[q-p + m 2 N ± WatM^] 



(IV.13) 



V, = 2 



+ 2 



[2(Q 2 - 9 • p) 2 + M 2 R {3Q 2 - 4q ■ p)] + 2 



m 



A" 



C%C£q-p-CXCg(Q 2 -q-p) 
M R m N 



[2M R =f 2m N M R + Q 2 - q ■ p] + 2 



°5 °5 ~" 



(Q 2 -<?-p) 



(Q 2 -q-p) 



''N 



A m N CYCj 2 



(2M| ± 2m N M R + Q 2 -q- 



P 



q-p 



(IV. 14) 



These are the important structure functions for most of the kinematic region. There are two additional structure 
functions, whose contribution to the cross section is proportional to the square of the muon mass. 



Va = 



V\2 



M 2 R 



[{2q ■ p - Q 2 ){q ■ p + m 2 N ) - M 2 R {m 2 N ± 7TI N Mr)] 



M 2 R 



[{2q -p-Q 2 ){q-p + m 2 N ) - M 2 R {m 2 N T m N M R )] 



{cX) 2 {2q- P -Q 2 ) , (cl) 2 {q- P ) 2 , a\cx 



N 



^ N M\ 



■q-p 



N 



[q-p + m 2 N =p m N M R ] 



' 1 \(2q -p - Q 2 ){q -p + m 2 N =F 2m N M R ) + g • pM 2 R ] + 3 J [2g • p + m 2 ^ T 2™^^^ + M| + Q 2 ] 

m N M R 



tun Mr 



7 



( 5 ' Ml 



^(2q-p-Q 2 



m 



m 2 N M R 



, vM (Q 2 -?-p) 2 + g 2 M| 



+^4^5 2 _2 ^^/^ 



-(Af 2 + Q 2 -q-p) 



N 



[q-p + m 2 N ± toatMr] 



m N M R 



[{2q-p-Q 2 ){q-p + m 2 N ± 2m N M R ) + M 2 R q ■ p] 



C^C^(2q ■ p + m 2 N ± 2m N M R ) - ■ p( Q 2 + (M R ± m N " 2 



(IV.15) 



(c 3 y ) 2 + (c 3 A ) ; 

MI 



(CD 2 ^ (C 5 p 2 (Q 2 + M 2 ) ^ ,W 



c y c y 

q -p [q -p + m 2 N + M 2 R ] + -3 ^ q ■ p [q ■ p + {Mr T m N ) 2 + Q 2 ] 



m N M R 



''N 



m 2 N M 2 R 



N 



q-p[q-p + m 2 N =p m^Ma] H ^ • p [q ■ p + (Mr =F ^at) 2 ] 

m N M R 



+ 



m N M R 
(Cf) 2 



q-p[q-p+ (M R ± m w ) 2 ] + C^^- [2q ■ p + Q 2 + (M R ± m N ) 2 ] 



A" 



i? 



c/c, 



4 ° 6 Q 2 



m 



^6 
Mr 



(g-p-Q 2 ) 



[g -p + m 2 N ±m N M R ] 



(IV. 16) 



B. Resonance P 3 3(1232) 

The method of extracting the form factors from the 
hclicity amplitudes, described in the previous section is 
applicable to any resonance. The helicity amplitudes for 
the ^33(1232) resonance were calculated in a similar man- 
ner and we obtain the following results: 
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Comparing helicity amplitudes from Eqs. 
(|IV. 1 7|) . l|I V. 1 8|i . (|IV. 1 9|> with the available data Q 
allows us to extract the form factors 
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Form factors Co and 



are in agreement with 



those obtained in the magnetic dominance approxima- 
tion (which was used in all the previous papers on neu- 
trinoproduction). The agreement has 5% accuracy and 
at the same time the nonzero scalar helicity amplitude 
is described correctly. The fit of the proton helicity am- 
plitudes for the form factors from Eq. H1V.20|I is shown in 
FigQ] Electromagnetic neutron form factors and vector 
form factors for the neutrino-nucleon interactions can be 
calculated according to Eq. I|III.5(I . 

Axial form factors have already been discussed sev- 
eral times, the way to obtain them is illustrated in Ap- 
pendix H] Eq. HA.4|) . The result is 
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A practical aspect with this resonance concerns the 
cross section of the tau neutrino interactions, which is 
discussed in Section [V] 
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FIG. 4: Helicity amplitudes for the p33(1232) resonance, cal- 
culated with the form factors from Ea. l|IV.20fl . Data are from 
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C. Resonance Pn (1440) 



At nonzero Q 2 data on helicity amplitudes for the 
P\i(1440) are available only for proton. Unlike the other 
resonances, the accuracy of data in this case is low and 
numerical values, provided by different groups differ sig- 
nificantly, as illustrated in Fig|S] In this situation we use 
for our fit only the recent data |25l l3fj| . 
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For spin-1/2 resonances the parametrization of the 
weak vertex for the resonance production is simpler 
than for the spin-3/2 resonances and is similar to the 
parametrization for quasi-elastic scattering. 

The matrix elements of the Pn resonance production 
can be written as follows: 
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where we use the standard notation <j vp = | [7", j p ] and 
the kinematic factors are scaled by /1 = mjf + M R in 
order to make them dimensionless. 

To extract the vector form factors we use the same pro- 
cedure as before and calculate the helicity amplitudes for 
the virtual photoproduction process. Since the resonance 
has spin 1/2, only the A x / 2 and S1/2 amplitudes occur: 
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FIG. 5: Helicity amplitudes for the Pn(1440) resonance, cal- 
culated with the form factors from Ea. llIV.240 , For A 1 / 2 the 
ata are from: |2ql (unshaded circles), |25{| (unshaded pen- 
tagon s'). |3fl| (full circles); for 51/2: I2a | ( unshaded up trian- 
gles), |25(| (unshaded down triangles), |3Cl| (full triangles) 



Matching the data against Eqs. (|IV.22|) . (|IV.23jl al- 
lows us to parametrize the proton electromagnetic form 
factors as follows: 
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The difference among the reported helicity amplitudes 
are larger than the estimated contribution of the isoscalar 
part of the electromagnetic current. For this reason we 
shall assume that the isoscalar contribution is negligibly 
small and use the relation A 
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contribution in the neutrinoproduction is now given as 
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The differential cross section is expressed again with 

The hadronic structure func- 



the general formula IjlV.l 
tions are calculated explicitly to be: 
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The helicity amplitudes 
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are used to extract the electromagnetic form factors. 

As in the case of Pn(1440) resonance, we choose here 
to fit only proton data on helicity amplitudes and neglect 
the isoscalar contribution to the electromagnetic current 
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where the upper sign corresponds to P\\ and the lower 
sign to S\i resonance. 

As it is shown in AppendixEI PCAC allows us to relate 
the two axial form factors and fix their values at Q 2 = 0: 

A {P) = i,M R + m N )m N a (p) a(p )(0) = _ 051 
A Q 2 +ml 1 1 w 
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FIG. 6: Helicity amplitudes for the Pn(1440) resonance, cal- 
culated with the form factors from Eq.I IIV.33l . For A1/2 data 
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The Q 2 dependence of the form factors cannot be deter- 
mined by general theoretical consideration and will have 
to be extracted from the experimental data. We again 
suppose, that the form factors are modified dipoles 
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D. Resonance 5*11(1535) 
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For the Sn the amplitude of resonance production is 
similar to that for P\\ with the 75 matrix exchanged be- 
tween the vector and the axial parts 



{s 11 \r\N) = u(p l ) 



n V n A ' 

92 ■ up A v 9z v 

zcr p g p75 - g{ 7 q 



u{p). (IV.31) 



We obtain the form factors 
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Notice here, that g\ falls down slower than dipole (at 
least for Q 2 < 3.5 GeV 2 ), supplying the most prominent 
contribution among the three isospin-1/2 resonances dis- 
cussed in this paper. This means experimentally, that 
the relative role of the second resonance region increases 
with increasing of Q 2 . Values of Q 2 = 1 — 2 GeV 2 
are accessible (and are not suppressed kinematically) for 
E v > 1.5 — 2 GeV. For these energies the 511(1535) and 
Di3(1520) resonances are observable in the differential 
cross section. 
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The axial form factors are determined from PCAC as 
described in the Appendix 1X1 
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The Q 2 dependence is again taken 
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We adopt this functional form, but one must keep open 
the possibility that it may change when experimental re- 
sults become available. 



V. CROSS SECTION FOR THE TAU 
NEUTRINOS 

Before we describe numerical results for the cross sec- 
tions in the second resonance region, we shortly dis- 
cuss the cross section for r-neutrinos and the accuracy 
achieved in different calculations. 

Recently we calculated the cross section of the reso- 
nance production 0] by taking into account the effects 
from the nonzero mass of the outgoing leptons. They 
generally decrease the cross section at small Q 2 . For the 
muon the effect is noticeable in the Q 2 — dependence of 
the differential cross section, but is very small in the in- 
tegrated cross section. 

It was shown, that the cross section is decreased at 
small Q 2 via 1) reduction of the available phase space; 
2) nonzero contributions of the W4 and W5 structure 
functions. Following [10|, we'll refer to the latter effect 
as "dynamic correction". To date, several Monte Carlo 
Neutrino Simulators use the Rein-Sehgal model @ of the 
resonance production as an input. In this model the lep- 
ton mass is not included. Thus, in Monte Carlo simula- 
tions the phase space is restricted simply by kinematics, 
but they do not take into account effects from W4 and W5 
structure function. Some calculations are also available, 
where the partonic values for the structure functions 

W 4 = (a), W 5 = W 2 ■ (q ■ p)/Q 2 (6) (V.l) 

are included. We compare these two approximations with 
our full calculations for the integrated and differential 
cross sections. Fig. [7| shows the integrated cross section 
for the reaction v T p — > fj," A ++ . 

One can easily see that taking the partonic limit for 
the structure functions is a good approximation in this 
case. Ignoring the structure functions, however, leads to 
a 100% overestimate of the cross section which is inac- 
curate. In both cases the difference comes mainly from 
the "low" (close to the threshold) Q 2 region, as it is il- 
lustrated in Fig. |H| where the differential cross section 
for the tau neutrino energy E v = 5 GeV is shown. One 
easily sees, that the discrepancy at low Q 2 reaches 30% 
for the partonic structure functions and more than 100% 
when W4 and W5 are ignored. 
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FIG. 7: The integrated cross section a(E) as function of the 
neutrino energy. The dashed curve includes phase space cor- 
rections only; the dotted curve includes phase space correc- 
tions and structure functions from the parton model; the solid 
curve includes phase space and structure functions calculated 
in our model. 
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v r p — » yT A ++ at E y = 5 GeV. The meaning of the curves is 
the same as in Fig. [7| 



FIG. 



VI. CROSS SECTIONS IN THE SECOND 
RESONANCE REGION 

We finally return to the second resonance region and 
use the isovector form factors to calculate the cross 
section for neutrinoproduction of the resonances. We 
specialize to the final channels vn —> R — > pir° and 
vn — ► R — > nir + , where both I = 3/2 and / = 1/2 
resonances contribute. The data that we use is from 
the ANL mEl, SKAT IH and BNL experiments. 
The ANL and BNL experiments were carried on Hydro- 
gen and Deuterium targets, while the SKAT experiment 
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FIG. 10: Differential cross section da/dW for the one-pion 
neutrinoproduction for BEBC experiment 



used freon (CF^Br). The experiments use different neu- 
trino spectra, there is, however, an overlap region for 
E v < 2.0 GeV. The data show that the BNL points are 
consistently higher that those of the other two experi- 
ments (see figure 4a, b in ref. ^3)- This is also evident 
in earlier compilations of the data. For instance, Sakuda 
[3l| used the BNL data and his cross sections are larger 
that those of Paschos et al.|3^ where ANL and SKAT 
data were used. A recent article 0] uses data from a 
single experiment ^(| , where the differences between the 
experimental results is not evident. The error bars in 
these early experiments are rather large and it should be 
the task of the next experiments to improve them and 
settle the issue. 

The differential cross section da/dW was re por ted in 
several experiments (see figures 4 in 16j, 1 in |17|. 4 in 
[l8| . 7 in H3). We plot the differential cross section 
da/dW in figure for incoming neutrino energies E v 
\ . 2 and 3 GeV. We note, that the second resonance peak 
grows faster than the first one with neutrino energy and 
becomes more pronounced at E v = 3 GeV. 

In Fig^|we plot our theoretical curves together with 
the experimental data from the BEBC experiment [2(j 
for (E u ) = 54 GeV. The theoretical curve clearly shows 
two peaks with comparable areas under the peaks. The 
experimental points are of the same order of magnitude 
and follow general trends of our curves, but are not ac- 
curate enough to resolve two resonant peaks. 

The spectra of the invariant mass are also plotted in 
figure 4 in Ref. [3 up to W = 2.0 GeV, but there is no 
evident peak at 1.4 < W < 2.0 GeV, in spite of the fact 
that the number of events is large. This result together 
with the fact that the integrated cross sections for pir° 



and 



arc within errors comparable suggest that the 



7=1/2 and 7 = 3/2 amplitudes are comparable. 

We study next the integrated cross sections for the final 
states ^~p7r° and n~mr + as functions of the neutrino 
energy. The solid curves in Fig.llllshow the theoretically 
calculated cross sections with the cut W < 2.0 GeV and 



the dashed curve with the cut W < 1.6 GeV. For pir° 
the solid curve goes through most of the experimental 
points except for those of the BNL experiment, which are 
consistently higher than those of the other experiments. 
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For the wr + channel our curve is a little lower than the 
experimental points. This means that there are contri- 
butions from higher resonances or additional axial form 
factors. Another possibility is to add a smooth back- 
ground which grows with energy. An incoherent isospin- 
1/2 background of approximately 5 • 10~ A0 (E U /I GeV — 
0.28) 1 / 4 cm 2 would be sufficient to fit the data, as it 
is shown by a double-dashed curve. By isospin con- 
servation, the background for the p7r° channel is deter- 
mined to be half as big. Including this background, which 
may originate from various sources, produces the double- 
dashed curves in Fig. 1111 Since experimental points are 
not consistent with each other, it is premature for us to 
speculate on the additional terms. 

For the P33(1232) the elasticity is high, but for the 
other resonances is ~ 0.5, which implies substantial de- 
cays to multipion final states. We computed in our for- 
malism the integrated cross section for multipion pro- 
duction. The results are shown in Fig. ^] with two cuts 
W < 1.6 GeV and W < 2 GeV. The experimental points 
are from Ref. |l6|. 
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We obtain values for two axial form factors by applying 
PC AC (see Appendix [XJ whenever the decay width and 
elasticity is known. For the spin-3/2 resonances there is 
still freedom for two additional axial form factors whose 
contribution may be important. This should be tested in 
the experiments. 

We present differential and integrated cross sections in 
Sections E| and |VI| For the P 33 (1232) we point out, that 
the structure functions W4 and W5 are important for 
experiments with r— leptons because they modify the Q 2 
dependence and influence the integrated cross section. 

The second resonance region has a noticeable peak in 
da/dW (Figjnj, which grows as E v increases from 1 to 
3 GeV. The integrated cross section for the I = 1/2 
channel also grows with energy of the beam and may 
require stronger contribution from the resonance region 
and a non-resonant background fFig llll) . 

Our results are important for the new oscillation ex- 
periments. In addition to the production of the reso- 
nances and the decays to one pion and a nucleon, there 
are also decays to two and more pions. Multipion de- 
cays contribute to the integrated cross section with the 
cut W < 2 GeV at the level of (2 - 3) • 1CT 39 cm 2 for 
E v > 4 GeV. Thus our results are useful in understand- 
ing the second resonance region and may point the way 
how the resonances sum up to merge at higher Q 2 into 
the DIS region. 
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VII. CONCLUSIONS 

We described in this article a general formalism for 
analysing the excitation of resonances by neutrinos. We 
adopt a notation for the cross section very similar to 
that of DIS by introducing structure functions. We give 
explicit formulas for the structure functions in terms of 
form factors. The form factors describe the structure of 
the transition amplitudes from nucleons to resonances. 
The vector components appear in electro- and neutrino- 
production. We use recent data on hclicity amplitudes 
from JLAB and the Mainz accelerator to determine the 
form factors including Q 2 dependences. We found out, 
that several of them fall slower than the dipole form fac- 
tor, at least for Q 2 < 2 — 3 GeV 2 . The accuracy of these 
results is illustrated in Figures |3 El LH ED 



APPENDIX A: DECAYS OF THE RESONANCES 
AND PCAC 

1. P 33 (1232) 

For P33(1232) isospin invariance defines the following 
effective Lagrangian for the ANir interactions: 
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The total width of the irN decay of each A++, A+, A 
or A - is calculated in a straight-forward way 
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where the pion momentum for the on-mass-shell reso- 
nance is 



1 



Pit 



2M h 



(Mj| — m 2 N — m 2 ) 2 — Am 2 N m\. 



For the experimental value Ta = 0.114 GeV, we obtain 
9a = 15-3 GeV -1 . The resonance width IjA.ll) is propor- 
tional to the third power of the pion momentum, so for 
the running resonance width we use 
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According to PCAC 
(R+\d^(0)\n) =-imlU 



where f„ = 0.97771^. 

For P 33 (1232) the relation (|A~2)l turns into 
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and we obtain in the limit m„ 
the two form factors 
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The denominator of the above formula is phenomenolog- 
ically extended as q 2 — m 2 ,. Making use of the relation 
(|A.3|I for Q 2 — > one also obtains = gAfn/Vi- Thus 
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For the A ++ the 7riVP vertex is V3 times bigger, so, 
strictly speaking, C A is also -\/3 times bigger. However, 
for historical reasons, the form factors are convention- 
ally defined for the vertex W + n — > R + and a factor y3 
appears in vertex for iy + n — > R ++ . 



The decay width to the irN is 

p 

The total width of the Pi3(1520) resonances is approx- 
imately 0.125 GeV and the elasticity is about 0.5. For 
this values we obtain go = 15.5 GeV -1 and the running 
width of the resonance is again proportional to the third 
power of the pion momentum. The PCAC relation turns 
into 
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which results in 
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with C A (D 13 ){Q 2 = 0) 



g D f n = -2.1. 



3. Pn(1440) 

For Pn(1440) the isospin-invariant Lagrangian is de- 
fined as 
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The decay width is 
rpn^TrJV 
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For the elasticity of 0.6 and the total width of 0.350 GeV 
we obtain gp = 10.9. 
The PCAC relation is 
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2. Di 3 (1520) 

For Z)i3(1520) the isospin-invariant Lagrangian of the 
P/i3(1520)V7r interactions is defined as: 
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At m w — > it leads to 
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(here the denominator is extended as before) and at 
Q 2 — > the coupling is 
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4. (1535) 



which at — > leads to 



For 5n(1535) the isospin-invariant Lagrangian is de- 
fined as 
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The decay width is 
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For the elasticity of 0.4 and the total width of 0.150 GeV 
we obtain gs = 1.12. 

The PCAC relation is again 
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a, Q n m N (M R -m N ) A 
5 3 (Sn) = g2 _ m 2 5i (Sii)- 



(here the denominator is extended as before) and at 
Q 2 — > the coupling is 



flr(^ii) 



2 ffs/x 



3 M R - m N 



= -0.21. 
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